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We calculate the Drell-Yan production cross sections and differential distributions in the transverse 
momentum and rapidity ol the J^'^ = 1 exotic hadrons 0(2170), X(4260) and Fi,(10890) at the 
hadron colliders LHC and the Tevatron. These hadrons are tetraquark (four-quark) candidates, with 
a hidden ss, cc and bb quark pair, respectively. In deriving the distributions and cross sections, we 
include the order as QCD corrections, resum the large logarithms in the small transverse momentum 
region in the impact-parameter formalism, and use the state of the art parton distribution functions. 
Taking into account the data on the production and decays of these vector hadrons from the e+e" 
experiments, we present the production rates for the processes pp(p) — >■ (^(2170)(— >■ 0(lO2O)7r'^7r^ — >■ 
K+K~-K+%~) + pp(p) X(4260)(^ J/V'tt+tt" fi+fi'-K+n') + and pp{p) -i> Y6(10890)(-> 
(T(15'), T(2S'), T(3S'))7r^7r^ — > fi'^ H^-k^tv^) + .... Their measurements at the hadron colliders will 
provide new experimental avenues to explore the underlying dynamics of these hadrons. 

PACS numbers: 14.40.Rt, 13.85.Ni 
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Exotic hadron spectroscopy now stands on firm foot- 
ing, thanks mainly to experiments during the last sev- 
eral years at the two e+e~ B factories, BaBar and Behe, 
which have reported an impressive number of such states 
in the mass region of the charmonia [l| . Most of these 
states defy a conventional cc charmonium interpretation, 
but their affinity to decay into the hidden charm states 
J/ip,ip' and into open charm states DD^*^ reveal that 
they have a cc component in their Fock space. Of partic- 
ular interest for us is the J^'^ — 1 state F(4260), dis- 
covered by BaBar Q in the initial state radiation (ISR) 
process e+e^ 7isRi^4260) — >■ 7isr J/i/'tt+tt", con- 
firmed later by CLEO Q and Belle 3, with the latter 
finding that two interfering Brcit-Wigner amplitudes to 
the J/tpTT^TT^ state describe the data better. Maiani et 
al. have interpreted F(4260) as the first orbital excita- 
tion of a diquark-antidiquark (tetraquark) state ([cs] [cs\). 
Particle Data Group (PDG) @ has assigned the name 
X(4260) for this resonance, which is what we also use. 

Evidence exists also for an ss state 1^^(2175) with 
the quantum numbers J^*^ = 1 , which was first ob- 
served by BaBar Q also in the ISR process e+e~ — > 
7isr/o(980)(/)(1020), where /q(980) is an 0++ scalar state, 
later confirmed by BESII Q and Belle 0. In [l^l, 
1^5(2175) is interpreted as a tetraquark [sg][sg] with one 
unit of relative angular momentum. This state is now 
called ^2170) by PDG Q, which we also use. Likewise, 
Belle [nl, [Til measured the state r6( 10890) in the pro- 
cess e+e- ^ n(10890) ^ (T(IS'), T(2S'), T(3S'))7r+7r-. 
The production cross sections and final state distribu- 
tions, in particular, the dipion invariant mass spectra, 
can be understood if Yb ^0890) is interpreted as a hidden 
hh tetraquark state |13l - [l5| . 

The aim of this Letter is to investigate the Drell-Yan 
production of the J^'~^ = 1~~ exotic hadrons at the Teva- 
tron and the LHC pp{p) ^ 7* V+..., with V being one 



of the states (/>(2170), A:(4260) or n(lG89G). The run- 
ning common thread is that all three arc candidates for 
the first orbital excitation of diquark-antidiquark states 
with a hidden ss, cc and bb quark content, respectively. 
Drell-Yan processes are theoretically better understood 
than the corresponding hadronic (prompt) production 
processes. Unfortunately, due to the very small leptonic 
branching ratios @, production of these exotic states in 
the traditional pair (£^ = e*, /i^) is not promising 
in the processes pp{p) — > 7* — s- V{—^ + .... 

We point out that the corresponding production 
cross sections are large enough to be measured at 
the LHC and the Tevatron, if, instead of the Icp- 
ton pair, one concentrates on the final states in 
which these exotic vector hadrons have been discov- 
ered in the e'^e~ annihilation experiments: 0(2170) — 
0(lO2O)/o(98O), X(4260) -> J/V'tt+tt", and n(10890) ^ 
(T(lS'),T(2S'),T(3S'))7r+7r-. The obvious advantage is 
that the essential input (branching ratios for the discov- 
ery channels times the respective leptonic widths) needed 
for estimating the cross sections, are all provided by the 
e+e~ experiments, yielding model-independent cross sec- 
tions irrespective of the nature of these states. On the 
other hand, these measurements are challenging due to 
the preponderance of the tt'^'tt" pairs from the underlying 
event in pp and pp collisions, and hence the combinato- 
rial background is expected to be formidable. However, 
we trust that, once the energy-momentum profile of the 
background tt+tt" pairs at the hadron colliders is well 
understood, the background can be effectively removed 
by appropriate cuts, enabling the experiments in carrying 
out significant measurements in this sector. 

The DY cross sections are based on the factorization 
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theorem (here X denotes a bunch of hadrons) 

a{pp/pp^ V + X) = j dxidx2y^Ja{xi)fb{x2) 

a.b 

y.a{a + h^V{p)+X), (1) 



where a, b denotes a generic parton inside a pro- 
ton/antiproton, V = 0(2170), X(4260), ^(10890) for the 
processes considered here with p = {p°,pt,P^) being 
the momentum 4- vector of the V, and fa{xi), fb{x2) arc 
the parton distribution functions (PDFs), which depend 
on the fractional momenta Xi{i — 1,2) (an additional 
scale-dependence is suppressed here). We shall adopt 
the MSTW (Martin-Stirhng-Thorne-Watt) PDFs ^ in 
our numerical calculations, and use another PDF set, the 
CTEQIO 17[, to estimate the uncertainties from this 



source. The process-dependent partonic cross sections 
a{a + — > + X) will be computed using the QCD 
perturbation theory. 

We recall that the leading order contribution comes 
from the sub-process gg — > 7* — >■ 1/ 
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with the color factor — 3. We include the leading or- 
der QCD (i.e., 0{as)) corrections, implemented following 
the pioneering papers [l^, [3] ■ This formalism is applied 
to calculate the differential distributions (Pa/dydp^, 
with the rapidity defined as y = ^ In j^z^ ■ The trans- 
verse momentum distribution at the tree level has the 
form 6{p'^). Perturbative QCD (gluon bremsstrahlung) 
generates a non-trivial p^-distribution. However, large 
logarithms of the type a" ln™(p^/py) arising from higher 
order QCD corrections spoil the perturbative expansion 
in the small transverse momentum region. These large 
logarithms must be resummed in order to improve the 
convergence of the perturbation theory. This is done in 
the Collins-Soper-Sterman (CSS) framework [20| where 
the resummation is carried out in the impact parame- 
ter space, yielding a simple form for the resummed px 
distribution 
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in which d^G^*^" / dydp^ reorganizes the singular terms in 
the ~> limit. Explicitly, this takes the form 
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with x\ = mv/\fsey, X2 = my / , and s is the 
square of the center-of-mass collision energy. The func- 
tion /{pt) in (pi) is introduced as a matching function for 
which wc use 1211 
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To estimate the uncertainty caused by the matching pro- 
cedure, we take Qmatch = (2±l)my, and this uncertainty 
will be included in the numerical estimates of the trans- 
verse momentum distributions. 

The Sudakov factor W{b, ^,p,xi,X2) is expressed as 



Wib, -,p,xi,X2) 
b 



dfl^ 
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\ni^A[g,{n)) + B{g,,{n)) 



and the coefficient functions A, B, Cqa/qa{xi/xi) and 
C qb / qb{x2 / X2) are expanded (in units of (as/Tr)^")). Some 



leading terms in these expansions are [2C 
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where the integration constants C'l , C2 in the Sudakov 
factor (not shown explicitly) and the renormalization 
scale /i in (jH) have been taken a.s Ci = fib = c = 2e'~^'^ 
and C2 = 1; where 7^ = 0.57722 is the Eulcr constant. 

The asymptotic term in ([3]) coincides with the pertur- 
bative results in the small pT region 



dydp^ 



dydp^ 
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so that in this region the resummed terms dominate. The 
factorization scale is chosen as /i = y^p'^ + my. 

As the large impact parameter 6 corresponds to a 
low momentum scale, a cutoff is introduced in the 
CSS formalism (20j . which replaces the parameter b by 
b* — b/y^l + fo2/6^j^^, with bt bounded from above 
by bmax- The non-perturbative effects to compensate 
this cutoff are incorporated into a phcnomenological 
function F^pib, my, xi,X2), and a commonly-adopted 
parametrization obtained by fitting the data on W and 
Z production [i^] is given by 
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where gi = 0.11 GeV^ 52 = 0.58 GeV^ 53 = 
-1.5 GeV"^ and Qo = 1.6 GeV for b^ax = 0.5GeV"^ It 
should be pointed out that the above value of Qo is not 
appropriate for (/)(2170), as in this case m^(^2no) < 2Qoj 
which would lead to an enhancement of the large 6-region 
instead of suppressing it, and therefore in our calculation 
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TABLE I. Masses, total and partial decay widths of the 
0(2170), X(4260) and ^(,(10890). Unless specified, all input 
values are taken from the PDG review 0| 

my (MeV) r (McV) V^^B (eV) 

(^(2170) 2175 ±15 61 ±18 2.5 ± 0.9 ^ 

X(4260) 4263t^ 108 ± 21 [4| 6.0li j^ [4] 

Yb(10890) 10888.4±|g [12] 30.7l^;? [12] 0.69l[i;2g ^ [12] 

S^^K+K- (48.9 ±0.5)% S^„(980)^.+ .- (50t^)%[23] 

^?.//V.^M+M- (5.93 ±0.06)% (2.48 ±0.05)% 

gT(2S)-.^+M- (1.93 ±0.17)% Bt(3S)^m+m^ (2.18 ±0.21)% 

^ Tee X B(</)(2170) -i- ?i(1020)/o(980)). 
Tee X B(X(4260) J/V'TT+TT-), Corresponding to Solution I. 

" Tee X B(yi,(10890) ^ T(15)7r+7r-) obtained from 
a = (2.78to Pt'. '^i> T(2S)7r+7r-, the cross section 
(4.82l:J;Jj}) pb gives TeeB = (1.20t°;*^) eV, while for 
Y(, — > T(35)7r"'"7r~ , the cross section (1.7lj]Q'3g) pb corresponds 
to TeeB = (0.42i;i;J|) eV. 



TABLE II. Cross sections (in units of pb) for the pro- 
cesses ppij)) 0(2170) (-> 0(lO2O)/o(98O) K+K--K+-K-), 
pp{p) — >■ X(4260)(— ;> J/V'tt+tt" — s> ^+^"7r+7r"), and pp(p) — >• 
F6(10890)(^ T(1S', 25', 35)^+71- At+zi^Tr+Tr-), at the 
Tevatron (^ = 1.96 TeV) and the LHC (x/s = 7 TeV and 14 
TeV), using the MSTW PDFs. A rapidity range {\y\ < 2.5) is 
assumed for the Tevatron experiments (CDF and DO) and for 
the LHC experiments (ATLAS and CMS); a rapidity range 
1.9 < y < 4.9 is used for the LHCb. 





0(2170) 


X(4260) 


Yb( 10890) 


Tevatron(|i/| < 2.5) 
LHC 7TeV {\y\ < 2.5) 
LHCb 7TeV (1.9 <y < 4.9) 

LHC 14TeV (|t/| < 2.5) 
LHCb 14TeV (1.9 < y < 4.9) 


r, O+0.9 

3.6111 

*.0-1.9 
o 7+1.9 
^. ' -1.6 


r, r,o+0.19 

".^^-0.09 

0.24«-? 
0.54t°;l^ 

o.3i«:?I 


n 0020+" '"'"^ 

U.UUZU_o.0005 

0040+° "°" 
0.0023l°:S 
0.0060l°;°°l° 
0.0033l°:°°li 



we use as input Qq ~ 1.0 GcV, which we adopt for the 
16(10890) and X(4260) eases as well. 

The eleetroniagnetic coupling constants gqqv are re- 
lated to the e+e~y coupling g^ev by gqqv = SqgeeV- 
The relevant e+e" experimental data which are used to 
derive these parameters arc collected in Table HI The 
entries for Tf.e{yb)B{Yb T{jiS)tt'^tt~) are obtained 
by using the relation ree(yb)B(Yb T{nS)TT~^TT~) — 
ry^my^cr(T(nS')7r+7r~)/(127r), with all three quantities 
on the r.h.s. taken from Belle [l^ . 

Having specified the formalism and the necessary in- 
puts, we present our numerical results. As the distribu- 
tions at the Tevatron and the LHC are rather similar, we 
show the figures only for the LHC. Rapidity distributions 
da/dy (in units of pb) for the three Drell-Yan processes 
at the LHC for ^ = 7 TeV are shown in Fig. [Ij (a) 
pp (0(2170) -> 0(lO2O)/o(98O) ^ K+K-tt+tt-) ± 
(b) pp -> (X(4260) ^ J/V'TT+TT- ^ fi+^rn+n-) + 
and (c) pp (n(10890) ^ T(1S', 25', 3S')7r+7r- ^ 



fi'^ fi~'K~^TT~) + ... (contributions from three intermediate 
states have been added). The normalized distributions 
are stable, though the indicated uncertainties in the nor- 
malization in Table HI] discussed below will also reflect in 
the rapidity distributions shown in this figure. The cor- 
responding transverse momentum distributions da/dpT 
(in units of pb/GeV) are shown in Fig. [2l which are ob- 
tained for the rapidity range \y\ < 2.5 (for ATLAS and 
CMS). The corresponding distributions in the rapidity 
range 1.9 < y < 4.9 (for the LHCb) are very similar, 
and hence not shown. The uncertainties caused by the 
matching functions are displayed. 

The integrated cross sections for the processes pp(p) — !• 
n(10890)(^ T(lS',2S',3S')7r+7r- ^ fi' tt+ tt- ) + 
PpIp) -> ^(4260)(^ J/tpTT+TT^ /i+^"7r+7r") + ..., and 
PpIp) 0(217O)(^ 0(lO2O)/o(98O) ^ K+K-tt+tt-) + 
...at the Tevatron (y^ = 1.96 TeV) and the LHC (for 
^/s = 7 TeV and 14 TeV) are presented in Table HH us- 
ing the MSTW PDFs [lal. The errors shown are from 
the parametric uncertainties in the PDFs and the vari- 
ous experimental inputs given in Table U which we have 
added in quadrature. We have also checked that our re- 
sults are modified only moderately if we use a different 
set of PDFs. For the CTEQIO PDFs [13, most changes 
amount to less than 30%, which are smaller than the 
uncertainties from the experimental input. We remark 
that the cross sections for CDF and DO (v^ = 1.96 TeV) 
and the LHCb (for = 7 TeV) are comparable, despite 
different center-of-mass energies, due to their different ra- 
pidity ranges, whereas the cross sections for the ATLAS 
and CMS detectors at the LHC are larger by typically 1.6 
(for 0(2170)), 1.7 (for A:(4260)) and 2.0 (for n(10890)), 
compared to the ones calculated for the CDF and DO at 
the Tevatron. Another remark concerns the collision en- 
ergy dependence. The cross sections at the LHC with 
= 14 TeV are enhanced by roughly 1.2 (for 0(2170)), 
1.3 (for X(4260)) and 1.5 (for n(10890)) compared to 
the corresponding results at y/s = 7TeV. 

To estimate the number of events, we assume an in- 
tegrated luminosity of 10 fb~^ at the Tevatron by the 
end of this year, and half that number at the LHC (for 
= 7 TeV) by the end of 2012. This yields 2.3 x 10^ 
events for the mode 0(2170) 0(lO2O)/o(98O) ^ 
K+K-TT+TT-, 2.3 X 10^ events for the mode X(4260) 
J/iJjTT+TT^ — > IJ.+ fJ.^ (and approximately the same num- 
ber for the Ar(4260) J/i/jtt+tt- — s- e+e~7r+7r~ mode), 
and only about 20 events for the mode Yt (10890) ^• 
(T(lS'),T(2S'),T(3S'))7r+7r- ^ fi+^rTT+TT- (and ap- 
proximately the same number of events for the 
n(10890) ^ (T(IS'), T(2S'), T(3S'))7r+7r- ^ e+e'TT+TT- 
mode). The corresponding numbers for the ATLAS and 
CMS [LHCb] are 1.8[1.1] x 10^, 2.0[1.2] x 10^, and 20 
[11], respectively. Hence, all these processes have mea- 
surable rates, given the luminosities at the Tevatron and 
the LHC, though the measurement of Yh(10890) in the 
Drell-Yan process may have to wait for higher luminosi- 
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y (a) y (b) y (c) 

FIG. 1. Rapidity distributions ^ (in units of pb) for the process (a) pp (^(2170) <?!>(1020)/o(980) -i> A'+A'"7r+7r") -I- 
(b) pp (X(4260) -> J/tpTv+Tv- ^i+fi-n+TT~) + and (c) pp ^ (^6(10890) T(1S, 25, 35)7r+7r" fi+ p,' n' ) + ... at 
LHC with = 7TeV using the MSTW PDFs. 




pAGeW) (a) PriGeW) (b) pAGeW) (c) 

FIG. 2. Transverse momentum distributions (in units of pb/GeV) for the process (a) pp — >• (0(2170) — >■ (j!)(1020)/o(980) — > 
K+K-TT+7V-) + (b) pp (X(4260) ^ J/V'tt+tt^ ^ /x+^-Tr+Tr^) -f and (c) pp (^(10890) T(15', 2S, 3S)7r+7r- ^ 
/x+^"7r+7r") + ... at the LHC ( = 7 TeV) with the rapidity cut \y\ < 2.5 using the MSTW PDFs. Uncertainties caused by 
the matching functions are displayed through 1/[1 -I- (pr/Qmatch)*] with Qmatch = (2 ± l)mv. 



ties and/or higher center-of-mass energy at the LHC. 

In summary, we have presented the Drell-Yan cross sec- 
tions and the corresponding differential distributions for 
the production of the J^'^ = 1 exotic vector hadrons 
(/)(2170), X(4260) and ^(10890) at the Tevatron and the 
LHC. The estimates given here arc model-independent 
due to the experimental input provided by the e~^e~ ex- 
periments. To unravel the dynamics underlying the ex- 
otic spectroscopy, one will have to undertake detailed dy- 
namical studies involving the final states. 

We acknowledge helpful discussions with Silja Brens- 
ing. Christian Hambrock and Satoshi Mishima. W. W. 
is supported by the Alexander- von-Humboldt Stiftung. 
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